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ABSTRACT 

We present the high resolution spectra of the young and compact planetary 
nebula NGC 6790 obtained with the echelle spectrograph at Bohyunsan Optical 
Astronomy Observatory and report the discovery of Raman scattered He II A 6545 
in this object. This line feature is formed in a thick neutral region surrounding 
the hot central star, where He IIA 1025 line photons are scattered inelastically 
by hydrogen atoms. A Monte Carlo technique is adopted to compute the line 
profiles with a simple geometric model, in which the neutral region is in the 
form of a cylindrical shell that is expanding from the central star. From our 
line profile analysis, the expansion velocity of the H I region lies in the range 
Vexp = 15 — 19 km s _1 . Less stringent constraints are put on the H I column 
density Nhi and covering factor C, where the total flux of Raman He IIA6545 is 
consistent with their product CNhi ~ 0.5 x 10 20 cm -2 . The Monte Carlo profiles 
from stationary emission models exhibit deficit in the wing parts. A much better 
fit is obtained when the He II emission region is assumed to take the form of a 
ring that slowly rotates with a rotation speed ~ 18 km s" 1 . Brief discussions are 
presented regarding the mass loss processes and future observations. 

Subject headings: planetary nebulae — planetary nebulae: individual NGC 6790 - 
radiative transfer — scattering — mass loss 
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1. Introduction 

Mass loss is an important process that mainly occurs in the late stage of stellar 
evolution. A star with a mass less than 8 M loses a significant amount of mass in the giant 
stage before becoming a planetary nebula with a hot white dwarf at its center. Considering 
the Chandrasekhar limit of 1.4 M , the mass loss process in the giant stage with enriched 
heavy elements should be important in the chemical evolution of the interstellar medium. 
In this regard, with a recent history of mass loss, young planetary nebulae are interesting 
objects to study the mass loss process. 

It is expected that around a young planetary nebula there may be a significant amount 
of neutral material that was lost in the previous stage of stellar evolution. In this case, the 
neutral region is exposed to the strong UV emission line source in the vicinity of the hot 
central star of the planetary nebula. Therefore, important information related with the 
mass loss process can be gathered from investigations of the scattering processes of the UV 
radiation originating from the center region. 

Taylor, Gussie & Pottasch (1990) performed H I 21 cm radio observations for a number 
of compact planetary nebulae (see also Altschuler et al. 1986, Gussie & Taylor 1995, 
Schneider et al. 1987). Their target selection was made on the basis of high radio brightness 
temperature, which is indicative of the nebular compactness. They searched an absorption 
trough that may be formed at the radial velocity of a compact planetary nebula when the 
neutral region blocks the background H I radio emission from our Galaxy. A number of 
compact young planetary nebulae including IC 5117 and NGC 6790 have been detected. 
Adopting an excitation temperature T HI = 100 K, the typical H I column density was 
determined to be of order N HI ~ 10 20 cm" 2 in these objects. 

Astrophysical Raman spectroscopy involving atomic hydrogen was initiated by Schmid 
(1989), who identified the mysterious broad emission bands occurring at 6825 A and 7088 
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A in many symbiotic stars (see also Nussbaumer, Schmid & Vogel 1989). He proposed 
that a hydrogen atom in the ground state is excited with the absorption of an incident far 
UV O VIA 1032 photon and de-excites into the 2s level with the re-emission of an optical 
photon at 6825 A. An analogous process for far UV O VIA1038 yields optical photons at 
7088 A. The large line width and prominent linear polarization exhibited by these scattered 
features strongly support his proposal (e.g. Harries & Howarth 1996). Observations made 
simultaneous in the UV and optical regions also confirm the Raman scattering nature 
(Espey et al. 1995). 

In the spectrum of the symbiotic star RR Telescopii, Van Groningen (1993) discovered 
Raman scattered He II features that are formed blueward of hydrogen Balmer emission 
lines. He II emission lines arising from transitions between n = 2k and n = 2 levels have 
wavelengths that are slightly shorter than hydrogen Lyman lines owing to the fact that 
He II ions are single electron atoms with a slightly larger two body reduced mass. The 
proximity to resonance is responsible for a large scattering cross section requiring the 
existence of a neutral region with Nhi ~ 10 20 cm -2 around a He II emission source. Raman 
scattered He II features are also reported in other symbiotic stars including He 2-106, 
HM Sagittae and V1016 Cygni (Lee, Kang & Byun 2001, Jung & Lee 2004b, Birriel 2004). 

Raman scattering of He II by atomic hydrogen also operates in young planetary 
nebulae. The first discovery was reported by Pequignot et al. (1997) in their spectroscopic 
analysis of the young planetary nebula NGC 7027. Subsequently, Groves et al. (2002) found 
the same Raman scattered He II features in the planetary nebula NGC 6302. Recently, 
Lee et al. (2006) reported that the compact planetary nebula IC 5117 also exhibits 
Raman scattered He II features blueward of Ha and H/5. In these objects, it appears that 
the central star is surrounded by a neutral region with a significant covering factor. In 
particular, Lee et al. (2006) discussed in detail the atomic physics of He II recombination 
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and Raman scattering processes. 

We present our high resolution spectra of the young and compact planetary nebula 
NGC 6790 and report our finding of the Raman scattered He IIA6545 feature in this object. 
Using the Ha image, Tylenda et al. (2003) measured the angular size of NGC 6790 to 
be 4" x 3". This size estimate of NGC 6790 is consistent with the HST image shown by 
Kwok, Su & Sahai (2003), who also identified two inner shells of similar orientations in 
NGC 6790. The distance to NGC 6790 is poorly known. Gathier et al. (1986) proposed 
that NGC 6790 is further than ~ 0.8 kpc based on their kinematic considerations. Adopting 
a statistical method Zhang (1995) suggested a distance of 5.7 kpc to NGC 6790. In their 
high resolution spectroscopy of NGC 6790, Aller, Hyung & Feibelman (1996) proposed a 
core mass of 0.6 M and an age of 6000 yr with the note that these values are dependent 
on the uncertain distance to NGC 6790. 

We perform Monte Carlo radiative transfer simulations in order to obtain the geometric 
and kinematic information of the neutral region. In section 2, we describe our observation 
and line fitting analyses and the following section presents our results of the Monte Carlo 
radiative transfer. In the final section, we discuss briefly our observation and mass loss 
processes of NGC 6790. 



2. Observation and Analysis 

2.1. Observation and Data 

We observed the young planetary nebula NGC 6790 on the night of 2008 May 
31 using the 1.8 m telescope at Bohyunsan Optical Astronomy Observatory (BOAO). 
The spectrograph that we used is the BOES (BOAO Echelle Spectrograph), which is a 
bench-mounted echelle system fed by optical fibers with various diameters. We used the 
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300 micron fiber, which yields the spectral resolution ~ 30, 000 with the field of view of 
3". The spectral coverage ranges 3600 A through 10,500 A. We obtained two spectra with 
exposure times of 600 s and 7200 s, respectively. A Th-Ar lamp was used for wavelength 
calibraions. For more detailed information on BOES, one is referred to Kim et al. (2007). 
Standard procedures using the IRAF packages were followed to reduce the spectra. 

In Fig. [H we show parts of our spectra around Ha and B./3. The vertical axis represents 
the relative flux density. We normalize the flux density using [N II]A6548, which is set to 
have a flux density peak of unity. The top panel of Fig. [T] is the spectrum around Ha with 
an exposure time of 600 s. We note strong forbidden emission lines of N II at 6548 A and 
6583 A. In this short exposure spectral image, the strongest Ha is unsaturated, allowing 
us to fit the Ha profile. The middle panel of Fig. [T] shows the Ha part of the spectrum 
with an exposure time of 7200 s. The strong Ha is saturated and we can discern very faint 
emission lines including He IIA6527. He IIA6527 arises from transitions between n = 14 and 
n = 5 levels. We also clearly notice that around [N II]A6548 there exists a broad bump-like 
feature. This feature is not an instrumental artifact because no such feature is present near 
[N II] A 6583, which is supposed to be 3 times stronger than [N II]A6548 (e.g. Osterbrock 
1987). We propose that this broad feature is Raman scattered He IIA6545. 

The bottom panel of Fig. (TJ we show our spectrum of NGC 6790 around Hf3 with 
the exposure time of 7200 s. If Raman scattered He II exists blueward of Ha, we may 
expect a similar feature blueward of H/5. Indeed, when Pequignot et al. (1997) reported 
the operation of He II Raman scattering in NGC 7027, they detected Raman scattered 
He IIA4850. In this object, Raman scattered He IIA4850 is not blended with other strong 
emission lines, which is in contrast with Raman He IIA6545 that is severely blended with 
[N II]A6548. In the bottom panel of FigdJ no broad feature around 4850 A is detected 
with a level of any significance. The quite strong and sharp emission feature at 4851 A 
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is an emission line totally irrelevant with Raman scattering. Aller et al. (1996) identified 
this emission line as a forbidden line from Fe II. Our spectrum is of insufficient quality to 
confirm the existence of Raman scattered He IIA 4850. However, this does not cast serious 
doubts of the Raman scattering nature of the 6545 feature, because Raman He IIA4850 
is always weaker than Raman He IIA6545. More discussion on this point is presented in 
section 3.2. 



2.2. Line Fitting Analysis 

Single Gaussian functions in the form /(A) = / ex P[ — (A — A C ) 2 /AA 2 ] are used to fit 
the permitted emission lines of Ha, He IIA6560, He IIA6527 and the two N II forbidden 
lines. The least chi square method is adopted to obtain the best fitting Gaussian functions. 
We use the atomic spectral data from the website of the National Institute of Standard and 
Technology (NIST), from which we note that each emission line in our spectra of NGC 6790 
appears systematically redward of atomic line center by an amount of 20.7 km s _1 . In Table 
1, we summarize the result of our profile analysis. The fitting parameters are quite similar 
to those found for IC 5117 by Lee et al. (2006). 

Fig. [2] illustrates our line fitting analysis of the emission lines in NGC 6790. The 
top panels show the result for Ha and He IIA6560. The short exposure data are used 
for the Ha emission line, which is excellently fitted by a single Gaussian function with a 
width A A = 0.54 A. He IIA6560 is also well fitted by a single Gaussian function with a 
considerably smaller width of AA = 0.48 A than that for Ha. 

The middle panels of Fig. [2] show our result for He IIA6527, which is significantly weak 
compared with He IIA 6560. He IIA6527 is strongly blended with another unidentified 
emission line. Because He IIA6560 is well-fitted by a single Gaussian, He IIA6527 should be 
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also fitted by a single Gaussian function, which is shown by the dotted line in panel (c). The 
long dashed line in panel (c) shows our Gaussian fit to the unidentified emission line. Groves 
et al.(2002) noted the existence of [N II]A6527 redward of He IIA6527 with the wavelength 
difference of 0.14 A in their spectrum of NGC 6302. However, the unidentified emission 
line in our spectrum of NGC 6790 can not be [N II] A6527, because it appears redward of 
He IIA6527 by 1 A. Furthermore, based on the NIST data, [N II] A6527 has the Einstein A 
coefficient A = 5.45 x lO^V 1 . Compared with [N II]A6548 having A = 9.19 x 10~ 4 , [N 
II]A6527 should be weaker than [N II]A6548 by a factor of 1700. Based on these atomic 
data, we plotted [N II] A6527 with a dot-dashed line in Fig.2. As is shown in the figure, [N 
II] A6527 is significantly weaker than He IIA6527, and hence can not affect the over all line 
fitting result. The 6528A feature is much stronger than [N II] A6527 and still remains to be 
identified. Panel (d) shows the composite profiles of the two single Gaussian functions in 
panel (c). From our profile analysis shown also in Table 1, we conclude that the flux ratio 
of He IIA 6527 and He IIA 6560 is 

^6527/^6560 = 4.1 X 10 " 2 . (1) 

The bottom panels of Fig. [2] show the detailed profiles of [N II] lines. It is interesting 
to note that [N II] A6548 exhibits a sharp absorption feature centered at 6548.60 A. The 
line center of [N II]A6548 appears at 6548.51 A, and the sharp absorption feature is 
excellently fitted by a single Gaussian function with a width of AA = 0.1 A and center 
at Ao = 6548.60 A. We find no such absorption feature in [N II] A6583, which should 
exhibit exactly the same profile with 3 times more flux (e.g. Osterbrock 1987). To our 
knowledge, no plausible metal transition is responsible for this sharp absorption. We also 
checked the telluric absorption lines without finding any strong candidate. In the spectrum 
of IC 5117 obtained with the 3.6 m Canada- France-Hawaii Telescope we find no similar 
absorption feature, for which [N II]A6548 exhibits exactly the same profile as [N II] A6583. 



- 9- 



We tentatively propose that this is attributed to Ha that is redshifted by an amount of 
v a bs ~ 800 km s _1 . However, in this work, we limit our attention to the Raman scattered 
He IIA6545 with no further discussion of this possibly interesting feature. 

Lee et al.(2001) performed a line profile analysis of Raman scattered He IIA6545 in a 
number of symbiotic stars. They subtracted one third of the flux near [N II] A6583 from 
the flux near [N II]A6548 to expose a broad Raman scattering line feature successfully. 
However, in view of the existence of the unidentified absorption feature in [N II]A6548 and 
more severe blending with [N II]A6548, we took another approach, in which the Raman 
scattered He IIA6545 feature is directly fitted from our Monte Carlo data. 

3. Monte Carlo Radiative Transfer 

3.1. Monte Carlo Procedure 

In this subsection, we describe the procedure of our Monte Carlo analysis of the Raman 
scattered He IIA6545. Many planetary nebulae exhibit nonspherical morphology, which 
may have its origin in the asymmetric mass loss processes. In the case of NGC 6790, the 
HST image obtained by Kwok et al. (2003) shows elongated shells around the central star. 
As a first approximation, we adopt a cylindrical shell model for neutral material, which is 
schematically illustrated in Fig. [31 A similar geometry was considered in the analysis of 
IC 5117 by Lee et al. (2006). 

In this cylindrical shell geometry, the hot UV source is located at the center and H I 
material is uniformly distributed inside the cylindrical shell with finite height and thickness. 
The same geometry was adopted by Lee et al. (2006). However, the essential difference 
is that we now consider the scattering region is expanding with the constant expansion 
velocity v exp . The cylindrical region is characterized by a uniform H I density uh, the 
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height H and the inner and outer radii R H and R H + AR, respectively. In this case, the H I 
column density of the cylindrical shell is given by N HI = n H AR. 

Since the shell is of uniform density, instead of the physical length / we measure the 
distance inside the shell in terms of the scattering optical depth r defined by 



where a to t is the sum of the cross sections for Rayleigh and Raman scattering. Since a tot 
is a sensitive function of a wavelength of the photon being considered, a given distance 
may correspond to different optical depths dependent on the wavelength. Therefore, once 
a photon is generated in the Monte Carlo simulation, we assume that the wavelength does 
not change as long as it is Rayleigh scattered. Considering that the scattering region is 
neutral, this assumption should be reasonable. 

The basic atomic physics of Raman scattering adopted in our Monte Carlo code is 
explained in detail by Jung & Lee (2004a). Due to the proximity of He II A 1025 to H I 
Ly(3 resonance, the scattering cross section increases steeply near Ly/3. Yoo, Bak & Lee 
(2002) showed that the branching ratio r b into Raman scattering increases approximately 
linearly with wavelength, which is given by 



where <jR am is the cross section for Raman scattering and A ^ is the Ly/5 center wavelength. 
Therefore, the Raman conversion into the optical region is quite sensitive to the incident 
wavelength, which in turn depends on the expansion velocity. 

From the energy conservation, a Raman scattered He II feature is characterized by its 
large width given by 



t = n H a tot l, 



(2) 



Tb — O Rami O tot 



0.1342 + 12.50(A - A L2//3 )/A L2//3 , 



(3) 




(4) 
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where Aj and A# am are wavelengths of the incident and Raman scattered radiation (e.g. 
Schmid 1989, Nussbaumer et al. 1989). In the case of Raman He II A6545, the profile width 
becomes about 6 times broader than He IIA1025, which endows a unique property that the 
profile is mainly determined from the relative motion between the emitter and the scatterer. 

In our Monte Carlo calculation, we also consider the re-entry of a photon emerging 
from the inner wall of the cylinder, for which we assume that this photon travel freely until 
it hits the inner wall on the opposite side. We consider a photon with a unit wavevector 
k supposed to travel a scattering optical depth r from the position = (xj, yi, Zi). If this 
photon emerges from the inner wall of the cylinder, we find the two points of intersection 
with the inner wall of the cylinder. This is accomplished by solving the quadratic for t p 

R 2 H = \(r l + T p k)-p\\ (5) 

for which we denote the two solutions by r pl and r p2 with r P 2 > r p \. Here, p is the unit 
vector pointing radially outward from the cylinder axis. The difference of the two solutions 
At p is given by 

Ar p = T p2 — T p i 

(1 - hi) ' W 

where k x , k y and k z are the components of k. By adding At p to the original photon path, 
we find the new scattering site in the other side of the shell. 

The incident He IIA1025 line flux and profile can be inferred from the case B 
recombination theory of single electron atoms provided by Storey & Hummer (1995). In 
Table 2, we show the expected He IIA1025 line flux relative to He IIA6560 and He IIA6527 
for electron number densities n e = 10 4 , 10 6 and 10 8 cm -3 and temperatures T e = 10 4 and 
2 x 10 4 K. We note that our observed flux ratio of He IIA 6527 and He IIA 6560 given 
in Eq. ([1]) is consistent with the nebular condition of n e ~ 10 6 cm -3 and T e = 10 4 K. 
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However, this choice is not unique and the range of He IIA1025 is already quite significant 
with the choice of parameters in Table 2. With this caveat in mind, we fix the electron 
number density n e = 10 6 cm -3 and T e = 10 4 K. Adopting these values of n e and T e , the 
recombination theory by Storey & Hummer (1995) gives -F1025 = 4.2F 6560 , which is used for 
our Monte Carlo calculations. 

The Monte Carlo simulation starts with a generation of He IIA1025 line photons having 
the same line profile with that of observed He IIA6560, and appropriately scaled using 
the recombination theory. As He IIA6560 is fitted by a single Gaussian with a width of 
AA = 0.48 A, we note that the line profile function fuy for He IIA1025 is given by 

fuvW = /l025 exp - [(A - Ai025)V A ^1025] ( 7 ) 

with AA1025 = 0.48 • 1025/6560 A = 0.075 A. Here, the peak value /1025 is appropriately 
adjusted to yield Fi 025 = 4.2F 6560 . 

We trace each individual He IIA1025 line photon until it escapes from the H I region. 
From Eq. it is noted that the profiles of the Raman scattered features are determined 
from the relative kinematics between the emission source and the H I region and almost 
independent of the observer's line of sight. Therefore, in this work, we collect all the 
photons irrespective of the final direction. 

3.2. Simulated Raman Profiles 

3.2.1. Spherical Emission Region 

In the work of Lee et al. (2006), the analysis of Raman scattered He IIA6545 was 
purely based on the atomic physics and focused on the exact location of line center. Their 
computation shows that the Raman scattered feature should be centered significantly 
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blueward of [N II]A6548. In Fig. [T], we note that the Raman He IIA6545 is completely 
blended with [N II]A6548, which implies that the neutral scattering region should be 
receding from the central UV source. 

In Fig. HI we show our Monte Carlo profiles for various expansion speeds v exp of the 
neutral scattering region with respect to the hot central star. In this figure, the height of the 
cylinder is taken to be infinite so that the covering factor of the scattering region is unity. 
The column density is fixed to Nhi = 1 x 10 20 cm -2 . The solid line shows our observed data 
and the other lines show our Monte Carlo profiles corresponding to various values of v exp . 
We can clearly notice the center shift of the Raman He IIA6545, which is highly enhanced 
due to the line broadening given in Eq. (j3J). The top panel shows the profiles for velocities 
Vexp < 40 km s" 1 . The bottom panel shows the profiles for velocities in the smaller range 
14 km s _1 < v exp < 22 km s -1 . From the figure, the plausible expansion velocity is around 
20 km s _1 , for which the peak wavelength resides inside the [NII]A6548 emission line. 

One interesting point to note from Fig. H] is that the strength of the Raman feature 
increases sharply as v exp increases despite the fact that the covering factor and Nhi are fixed. 
This is explained by the fact that the Raman scattering cross section sharply increases near 
Ha; due to Ly/3 resonance in the parent wavelength space. Therefore, a receding H I region 
yields more Raman scattered He IIA6545 photons than when the same region is stationary. 
This complicated dependence of the scattering cross section on wavelength also results in 
slightly asymmetric Raman profiles, which is barely noticeable in Fig. HI Therefore, the 
Raman conversion efficiency may be estimated accurately only after the kinematics of the 
scattering region with respect to the emission source is carefully determined. 

In the left panel of Fig. [5l we show the Raman profiles for various H I column 
densities ranging Nhi — 10 19 — 1.5 x 10 20 cm -2 with the fixed values of H/Rh = 2 and 
Vexp = 20 km s _1 . Within this range of Nhi, the overall strength is nearly proportional 
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to the H I column density, because the H I region is mostly optically thin with respect 
to Raman scattering of He IIA1025. This expansion speed is very similar to the value of 
16 km s _1 determined from Doppler shifted Na D absorption lines by Dinerstein, Sneden & 
Uglum (1995). 

The right panel of Fig. shows the Monte Carlo Raman profiles for various covering 
factors of the cylindrical shell. As is expected, the overall strength is also proportional to 
the covering factor. In both the panels of Fig. [51 we obtain qualitatively similar profiles. 
This implies that the Raman profile analysis severely suffers from the degeneracy problem 
involving the covering factor and H I column density. 

With this caveat in mind related with the degeneracy in Nhi and the covering factor, 
we show our best fit profile from the Monte Carlo calculations in Fig. [6J The model 
parameters are v exp = 19 km s _1 , N^i = 9 x 10 19 cm -2 and H/Rh = 1.7. As in IC 5117, 
the H I region significantly covers the hot central star in NGC 6790. However, in this figure 
we notice that the model profiles exhibit deficit both in the blue wing and red wing parts. If 
this deficit is real, then it implies that in the direction to the H I region the incident profile 
is broader than in the observer's line of sight. The next subsection discusses this point. 

Jung & Lee (2004b) developed a Monte Carlo code to compute the line profile of 
Raman scattered He II 4850 and analyzed their spectrum of the symbiotic star V1016 Cyg. 
Using the same code, we show in Fig. 7 the Monte Carlo profile for Raman scattered He II 
4850 by a long dashed line. The same column density and covering factor as in Fig. 6 were 
used in this calculation. In the figure, the solid line shows the BOES data with the exposure 
time of 7200 s. Our observational data are barely consistent with our interpretation of 
Raman scattering nature. The poor quality of the current observational data hinders a 
further serious quantitative analysis. A more fruitful analysis may be made only after 
observational data with a better quality are secured. 
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3.2.2. Ring-like Emission Region 

In this subsection, we perform line profile analyses in the case where the emission 
region takes the form of a ring that is rotating in the vicinity of the hot central star. In the 
previous section, it was assumed that the He II emission region is spherically symmetric 
and stationary. However, it is highly probable that the distribution of nebular material 
significantly deviates from spherical symmetry considering the non-spherical shape exhibited 
by most planetary nebulae (e.g. Corradi & Schwarz 1995). In this case, the emission region 
may plausibly possess an ordered motion component, which may also be associated with 
the nonspherical nebular morphology. Therefore, we may expect that ionized material is 
concentrated on the equatorial plane having some slow rotation velocity component. 

There exists little kinematic information available on the emission region very near the 
central star. No observational data of NGC 6790 are available in the archives of HUT and 
FUSE. In consideration of the absence of a unique kinematic model accounting for all the 
observed emission line profiles, we adopt a simple ring-like emission region, in which we 
investigate the line of sight effect on the profiles of the He II emission and Raman scattered 
lines. Depending on the line of sight of the observer, the rotation velocity component is 
reduced by the factor sinz, where i is the inclination angle of the ring. However, Eq. 
dictates that the Raman profile is determined by the velocity component of the emitter 
with respect to the scatterer and fairly insensitive to the line of sight. 

This proposition leads to an interesting interpretation of our profile fitting of Ha and 
He IIA6560 presented in the previous section. We may decompose the emission profiles into 
a bulk component and a random component. We further assume that the bulk component 
represents a slow rotation in the equatorial plane and that the random component is 
attributed to a thermal motion and a turbulent motion. For the sake of simplicity, we 
assume that He IIA6560 and Ha are formed in the same ring-like region that is in slow 
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rotation in the equatorial plane with the speed v bulk- 

A He ion being 4 times heavier than a hydrogen nucleus, the line width of He II 
due to the thermal motion is half of that for Ha if they are formed in the same region. 
However, if the emission region possesses some turbulent component, then overall random 
motion component for hydrogen is broader than that of He II by a factor less than 2. If we 
denote the electron temperature of NGC 6790 by T e = 10 4 T4 K, then the thermal velocity 
associated with Ha is given by 



Vth,H 



^ = 13 T 4 1/2 km s~\ 



where m p is the proton mass and ks is the Boltzmann constant (e.g. Rybicki & Lightman 
1979). Introducing v tur b for the turbulent velocity scale, we denote the random velocity 
components of Ha and He IIA6560 by v ran ,H and v ran ,He-, respectively, where 

Vran,H = V tu rb + V t h,H 
Vran,He = V turb + (v th>H / 2) . (9) 

Noting that there are three model parameters, namely i, Vbuik and v tur b for the two 
line widths, we also encounter a degeneracy problem. Hoping that future observations may 
provide independent constraints on some of these model parameters, we just pick out a set 
of values that yield a reasonable fit to our observed data. In the top panels of Fig. [HJ we 
show model line profiles for He IIA6560 and Ha from one such set consisting of 

sin i = 0.6, Vbuik = 18 km s _1 , v tur b = 14 km s _1 , 

V t h,H = 14 km S" 1 , V t h,He = 0.5v th ,H- (10) 

The thermal velocity v t h,u — 14 km s _1 is consistent with the electron temperature 
T e = 10 4 K, which is similar to that obtained by Aller et al. (1996) from their 
photoionization modeling. The overall fits to both Ha and He IIA6560 appear quite good. 
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The bulk velocity component is consistent with the size of the emission ring region of order 
1 AU if we interpret the bulk motion to be Keplerian. However, the bulk motion may not 
be related with the Keplerian motion but may be related with the rotation component of 
the central star, for which case the physical size of the emission region can be at best poorly 
constrained. 

Because the H I region is also concentrated on the equatorial plane, the full bulk velocity 
component should be considered without the inclination effect for far UV He IIA1025 that 
is incident on the H I region. In the bottom panel of Fig [HJ the dotted line shows the 
He IIA6560 profile that would be measured by a hypothetical observer in the equatorial 
plane. It is excellently fitted by a single Gaussian function with a width AA = 0.61 A, 
which is significantly larger than the observed value of AA = 0.48 A by a factor of 1.3. 
Hence, the emission profile for He A1025 incident on the neutral region should also be 
broadened by the same factor. 

In Fig. [9l we show our Monte Carlo result using the profile shown in the bottom 
panel of Fig. [S] and appropriately scaled to He IIA1025. The other model parameters 
are also adjusted for better fit and they are v exp = 15 km s _1 , Nhi = 9 x 10 19 cm -2 and 
H/Rh = 1.2. A much better fit is obtained than that considered in the previous section. 
However, it should also be pointed out that in constructing the profile in Fig. [9] more model 
parameters have been used than in the previous section and still the degenerate nature of 
the problem persists. 

The expansion velocity of the H I shell in Fig. [9] is only v exp = 15 km s -1 , which is 
significantly smaller than the value v exp = 19 km s _1 presented in Fig. 6. This notable 
discrepancy in expansion velocity is attributed to the scattering cross section that is sharply 
peaked around Ha. According to Jung & Lee (2004a), this leads to the center shift of a 
Raman scattered He II feature, which is dependent on the column density. The result shown 
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in Fig. [9] implies that the shape or the width of the incident profile also affects the location 
of the line center. A more quantitative investigation in a significantly large parameter space 
is left to the future work. 

4. Discussion 

H I Raman spectroscopy provides an accurate determination of the expansion velocity 
of the H I region, for the measurement of which H I 21 cm radio observation has been the 
unique tool so far. Our analysis shows that the expansion velocity lies between 15 km s _1 
and 19 km s _1 , which is consistent with the value of 16 km s _1 provided by Taylor et al. 
(1990). As was pointed out by Lee et al. (2006) the Raman spectroscopy allows one to 
determine the H I column density whereas the excitation temperature should be assumed 
before Nhi is deduced from H I 21 cm radio observation. According to Taylor et al. (1990), 
Nhi = 2.7 x 10 20 cm -2 assuming the excitation temperature Thi = 100 K. Our Raman 
profile analysis lends support to this excitation temperature. 

Our current data are of insufficient quality to lift the degeneracy of the covering factor 
and H I column density, and the overall strength of the Raman feature is determined from 
the product of the two quantities. However, our Monte Carlo calculations show that Raman 
profiles exhibit redward asymmetry due to enhanced scattering cross section toward Ha 
resonance. With better quality spectra that may be available from bigger telescopes, it is 
hoped that tighter constraints are obtained from more refined profile analyses. If Raman 
scattered He II 4850 blueward of H/3 can also be used, additional constraints can be put to 
break the degeneracy. 

Even though the distance to NGC 6790 is highly uncertain, we may assume that the 
distance is about 1 kpc for simple order of magnitude calculations. According to Tylenda 
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et al. (2003), the angular size of NGC 6790 is ~ 3". This gives a physical size of the H I 
region R ~ 5 x 10 16 cm. If the H I region is of a thin cylindrical shell with the height similar 
to its radius, the total number N tot of hydrogen atoms inside the shell is approximately 
given by N tot = 2irR 2 N H i ~ 6 x 10 54 . Here, in our order of magnitude estimate, we ignore 
the inclination effect, which will overestimate the total number of hydrogen atoms by the 
factor sin i. The H I mass of the neutral region is inferred to be M^i ~ 4 x 10~ 3 M Q . 

Furthermore, the expansion velocity of v exp ~ 15 km s" 1 and the physical size of 
R ~ 5 x 10 16 cm together imply the age of order of a thousand years for NGC 6790. It 
should be pointed out that these rough calculations are highly dependent on the assumed 
distance to NGC 6790 and still the physical size of the H I region is quite uncertain. 

The origin of sharp absorption feature that appeared in [N II]A6548 is quite uncertain. 
If this absorption feature is attributed to Ha, then it may imply the existence of clumpy 
components having a small covering factor with respect to the [N II] emission region and 
receding with a significant velocity of ~ 800 km s -1 . In some planetary nebulae including 
M2-9 and NGC 6543, it is known that fast collimated outflows exist around the central star 
with a velocity of order 1000 km s _1 (Balick 1989, Gruendl, Chu & Guerrero 2004, Prinja 
et al. 2007). Ueta, Fong & Meixner (2001) presented near IR imaging observations of 
AFGL 618 and reported their findings of molecular bullet-like features moving faster than 
200 km s -1 . However, it still remains a mystery whether a clumpy bullet-like object can be 
ejected with so large a velocity from the center region. 

It should be pointed out that a ring-like emission model may not be a unique choice 
for the observed profiles of He IIA6560 and Ha. Many kinematical models involving 
jet-like outflows or radial infall and/or outflows may also yield similarly well-fitting profiles. 
Therefore without convincing support from other studies such as imaging observations using 
interferometry or hydrodynamical computations, it appears to be too early to conclude 
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about the kinematics of the He II emission region. 

A ring-like emission region and H I region concentrated in the equatorial region 
may provide interesting opportunities for spectropolarimetry. In symbiotic stars, Raman 
scattered O VIAA6825, 7088 are known to exhibit strong linear polarization (e.g. Harries 
& Howarth 1996, Schmid 1998). The polarization structure may be closely related with 
the accretion and mass loss processes that deviate from spherical symmetry (e.g. Lee & 
Park 1999, Lee & Kang 2007, Ikeda et al. 2004). Because Raman scattered features consist 
of purely scattered photons, they make ideal targets for linear spectropolarimetry. Future 
spectropolarimetric studies may provide more interesting information regarding the mass 
loss processes in AGB stars and planetary nebulae. 
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Table 1. Single Gaussian Fit Parameters of Emission Lines 



Line 


Ao (A) 


fo 


AA (A) 


Ha 6563 


6563.23 


34.8 


0.54 


He II A 6560 


6560.58 


0.072 


0.48 


He II A 6527 


6527.49 


0.00295 


0.48 


[N II] A 6548 


6548.51 


0.897 


0.47 


[N II] A 6583 


6583.90 


2.73 


0.48 



x 30 

a; 20 

^> 

iS 10 

0) 



r NGC 6790 

- [N II] , 

6548 


iii 

\ Ha ^ 
6563 ; 

[N M] i 
6583 : 

I _. _ A - 


- i i i i i i i i i i - 


6540 6560 


6580 

i i i i 





0.1 


X 


2 




0) 




> 


0.05 






a 

Sh 





He II 

6527 



[NM] 
6548 



Raman He II 
6545^ 

_L 



He 
6560 




Ha 
6563 



[N II] 

6583 




6540 



6560 



6580 



0.03 



0.02 



£ 0.01 








4850 



4860 



4870 



x(A) 



Fig. 1. — High resolution spectra of NGC 6790 obtained with BOES. The top panel is a short 
exposure spectrum with exposure time of 600 s, and the exposure time for the middle and 
bottom panels is 7200 s. The relative flux density is normalized such that [N II] A6548 has the 
flux density peak of unity. In the top panel, the strongest emission line Ha is unsaturated. 
We also clearly see [N II] lines at 6548 A and 6583 A. In the middle panel, we note that 
around [N II] A6548 there exists a broad wing feature. No similar feature is present around 3 
times stronger [N II] A6583, which means that the broad wing feature around [N II] A6548 is 
not associated with [N II] nor is an instrumental artifacts. The bottom panel shows the H/3 
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Fig. 2. — Gaussian line fitting analysis. The solid lines show the observational data and the 
dotted lines show our Gaussian fits. The same flux normalization as in Fig.l is used. The 
top panels show the results for Ha and He IIA6560. The middle panels show the line fitting 
result for He IIA6527 and a nearby unidentified emission line. He IIA6527 is fitted by a single 
Gaussian function, which is shown by the dotted line in panel (c). The long dashed line in 
panel (c) shows the unidentified emission line. Using the atomic data provided by NIST, we 
show the line contribution from [N II]A6527 by a dot-dashed line. In panel (d) we show the 
composite profile from the three single Gaussians shown in panel (c). The bottom panels 
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Fig. 3. — A schematic diagram of the Raman scattering geometry adopted in this work. The 
hot star and He II emission region are located at the center. Surrounding the UV emission 
region, the H I scattering region takes the form of a cylindrical shell with the inner radius 
Rh, the outer radius Rh + Ai? and the height H . In this work, the cylindrical shell is 
assumed to expand with the speed v exp . Hydrogen atoms are distributed uniformly with a 
number density nn inside the cylindrical shell. 
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Fig. 4. — Line profiles of Raman scattered He IIA6545 from our Monte Carlo simulations 
for various expansion speeds. The covering factor is fixed to be unity and Nhi = 10 20 cm -2 . 
Due to the inelasticity of Raman scattering or Eq.(jH), the location of line center is fairly 
sensitive to v exp . The top panel shows the profiles for velocities in the range v exp < 40 km s _1 
in an interval of 10 km s _1 . The bottom panel shows the profiles for velocities in the range 
14 km s _1 < v exp < 22 km s _1 . It is notable that the expansion velocity v exp affects both 
the location of line center and the total Raman flux. 
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Fig. 5. — Monte Carlo line profiles of Raman scattered He IIA6545 for various N HI and 
covering factors. The left panel shows the Monte Carlo profiles for various Nhi with the 
covering factor fixed to be H/R H = 2. The right panel shows the simulated profiles for 
various covering factors with fixed N HI = 10 20 cm -2 . 
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Fig. 6. — Our best fit Monte Carlo profile of Raman scattered He IIA6545 from a stationary- 
emission region surrounded by a cylindrical shell. The dotted line is the Monte Carlo line 
profile and the solid line is the observed data. The adopted parameters are v exp = 19 km s _1 , 



H/R H = 1.7, N m = 9 x 10 19 cm" 2 . 
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Fig. 7. — BOES data around H/3(solid line) and the Monte Carlo profile of Raman scattered 
He IIA4850 (long dashed line). The same column density and covering factor as in Fig. 6 were 
used in the Monte Carlo calculation. The observational data are barely consistent with the 
Monte Carlo result. 
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Fig. 8. — Line profiles of He IIA6560 and Ha from a ring-like emission region. The axis of 
the ring makes an angle % with the line of sight, where we take sini = 0.6 as an example. 
The upper panels show line profiles of He IIA6560 and Ha viewed from the observer's line 
of sight. The fitting parameters are Vbuik = 18 km s" 1 ,^ = 14 km s" 1 and v t h,H = 
14 km s _1 , v t h,He = 0.5v t h,H- See the text of the definitions of these velocities. The solid lines 
represent the BOES data and the dotted lines are model profiles. The lower panel shows the 
observed He IIA6560 profile (solid line) and the model profile that would be observed in the 
equatorial direction. The dotted model profile is excellently fitted by a single Gaussian with 

n wirltii A X = fi1 A 
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Fig. 9. — A Monte Carlo best fit profile (dotted line) of Raman scattered He IIA6545 from a 
ring-like emission region considered in Fig. 7. The adopted model parameters are sini = 0.6, 
Vbuik = 18 km s , 

v turb — 14 km s 1 , Vth h — 14 km s 1 . Refer the text for the definitions 
of these parameters. This profile provides a much better fit than that shown in Fig. 6. It is 
noted that the expansion velocity of the H I shell is v exp = 15 km s _1 , which is significantly 
smaller than that considered in Fig. 6. 
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Table 2. He II Recombination Data by Storey & Hummer (1995) 



Line Ratio T e = 10 4 K T e = 2 X 10 4 K 



n e = 10 4 cm 3 

^1025/^6560 3.600 4.519 

^6527/^6560 3.952 X 10~ 2 4.085 X 10~ 2 
n e = 10 6 cm- 3 

^1025/-F6560 3.804 4.676 

^6527/^6560 4.098 X 10~ 2 4.152 x 10~ 2 
n e = 10 8 cm -3 

^1025/^6560 4.439 5.181 

^6527/^6560 4.942 X lO" 2 4.614 x 10- 2 



